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We analyze the Belle data [K. F. Chen et al. (Belle Collaboration), Phys. Rev. Lett. 100, 
112001 (2008); I. Adachi et al. (Belle Collaboration), arXiv:0808.2445] on the processes e+e~ -> 
Y(lS) 7r + 7r~, T(2S) 7r + 7r~ near the peak of the T(5S) resonance, which are found to be anoma- 
lously large in rates compared to similar dipion transitions between the lower T resonances. Assum- 
ing these final states arise from the production and decays of the J PC — 1 state Y b ( 10890), which 
we interpret as a bound (diquark-antidiquark) tetraquark state [6g][6q], a dynamical model for the decays 
Y(, —5- T(1S) 7r + 7r~, T(2S) 7r + 7r~ is presented. Depending on the phase space, these decays receive 
significant contributions from the scalar ++ states, /o(600) and /o(980), and from the 2 ++ qq-meson 
/2(1270). Our model provides excellent fits for the decay distributions, supporting Y b as a tetraquark state. 



The observation of the T(IS') vr + vr~ and T(2S) vr + vr~ 
states near the T(5S I ) resonance peak at ^/s = 10.87 
GeV at the KEKB e + e~ collider by the Belle collabo- 
ration d|] has received a lot of theoretical attention [Q]. 
The two puzzling features of these data are that, if in- 
terpreted in terms of the processes e + e _ —> T(5S) —> 
T(liS) vr + vr~, T(2S) vr + vr~, the rates are anomalously 
larger (by more than two orders of magnitude) than the 
expectations from scaling the comparable T(4S) decays 
to the T(5/S), and the shapes of the distributions in the 
dipion invariant mass m W7r and the cosine of the helic- 
ity angle, cos#, where 6 is the angle between the ir~ 
and T(5iS) in the dipion rest frame, are not described 
by the models J3D based on the QCD multipole expan- 
sion Q4J, [5|] - a feature also at variance with similar dipion 
transitions between lower T resonances. A critical obser- 
vation towards understanding these features is that the fi- 
nal states in question are produced not from the decays 
of T(55), but from the process e+e~ -)• Y b ( 10890) -)• 
T(15) vr+7r-,T(25) vr+vr", with Y b a 1~ state, hav- 
ing a total decay width T(Y b ) = 55 ± 9 MeV Jg]. In 
a closely related recent paper 0], we have analyzed the 
BaBar data [8] obtained at the SLAC B factory during 
an energy scan of the e + e" — > bb cross section in the 
range of the center of mass energy ^/s = 10.54 to 11.20 
GeV, observing that the BaBar data on the i? b -scan are 
consistent with the presence of an additional bb state Yn, q ] 
with a mass of 10.90 GeV and a width of about 30 MeV, 
apart from the T(5S) and T(6S) resonances. Identify- 
ing the J PC = 1 state Yj^j (10900) seen in the en- 
ergy scan of the e + e~ — > bb cross section by BaBar [8J] 
with the state Y b (10890) seen by Belle O," 

we present 

a dynamical model based on the tetraquark interpreta- 



tion of Y b ( 10890) and show that it is in excellent agree- 
ment with the measured distributions in the decays Y b — > 

T(15) TT+TT-,T(2S) 7T+7T-. 

In the tetraquark interpretation, Y[ bq ] is a J PC = I"" 
bound (diquark-antidiquark) state having the flavor con- 
tent Y[ bq ] = QQ = [bq][bq] (here q = u or q = d, 
and Q is a diquark) with the spin and orbital momen- 
tum quantum numbers: S Q = 0, Sq = 0, S q q = 
0) Lqq = 1 OJ. The first two quantum numbers are 
the diquark spin, antidiquark spin, respectively, and the 
last two denote the spin and the orbital angular quan- 
tum numbers of the tetraquarks, with the total spin being 
J = Sqq + Lqq = 1. Such spin-0 diquarks are called 
"good" diquarks itlfJl and an interpolating diquark oper- 
ator can be written as Q ia = e a07 (h^ 5 q] - gf c 7 5 ^ 7 ) 
(with qi = u,d for i = 1,2 and b c the charge con- 
jugate fc-quark field b c = — z6 T a" 2 75). The "good" di- 
quark Q ia is in the attractive anti-triplet (3) color chan- 
nel (with the color quantum numbers denoted by the Greek 
letters). There are two such J PC = 1~~ states, Yt bq ] = 
([bq]s=a[bq]s=o)p-wavc, with the mass eigenstates, called 
Yj* ji and Y[ b j t ] in @|, being orthogonal combinations of 
Yt bu i and Yn,s . Their mass difference is induced by isospin 
splitting m d — m u and a mixing angle and is estimated 
as AM(Y fc ) = (5.6 ± 2.8) MeV. In the following, we 
will not distinguish between the lighter and the heavier of 
these states and denote them by the common symbol Y b . 
The decays Y b — > T(1S) vr + vr~, T(2S) vr + vr~ are sub- 
dominant, but Zweig-allowed and involve essentially the 
quark rearrangements shown below. 

With the J PC of the Y b and T(raS) both 1~, the vr+vr" 
states in the decays Y b — > T(1S) vr + vr _ , T(2S) vr + vr~ 



2 



are allowed to have the ++ and 2++ quantum num- 
bers. There are only three low-lying states in the Parti- 
cle Data Group (PDG) [IIP which can contribute as inter- 
mediate states, namely the two ++ states, / (600) and 
/o(980), which, following IU2U13I1 . we take as the lowest 
tetraquark states, and the 2++ qq-meson state / 2 (1270), 
all of which decay dominantly into vrvr. For the decay 
Y b — > T(15) vr+7r~, all three states contribute. How- 
ever, kinematics allows only the / (600) in the decay 
Y b — > T(2S) vr+vr~. In addition, a non-resonant con- 
tribution with a significant D-wave fraction is required by 
the data on Y b -4 T(15) vr+vr-, T(2S") vr+vr^ . The dy- 
namical model described below encodes all these features. 
We start by showing the relevant diagrams for the decays 

Y b (q)^T( P )+7T + (k 1 )+TT-(k 2 ). 



(1) 



fo(i 



The initial state represents the tetraquark states Y b = 
[bq][bq], and T stands for T(15) and T(25). Both di- 
agrams involve the creation of a qq pair from the vac- 
uum, with diagram a resulting into the (non-resonant) fi- 
nal states T(15) vr+vr~ and T(2S) vr+vr~, and diagram 
b leading to the final states T(1S) (/„ (600), /„ (980)) 
and T(2S) / (600), with the implied subsequent decays 
(fo (600), f a (980)) -4 vr+vr-. The 2++ intermediate state 
/ 2 (1270) contributing to the decay Y b — > T(1S) vr+vr - is 
depicted below. 



(2) 



£.(1270) 



Writing the Lorentz -invariant amplitudes as 

M=e Y (q)el(p) V M?" (p, k u k 2 ) , (3) 



i=a,b 1 c 



where £^{q) and £^{p) are the polarization vectors of the 
Y b and T(nS), respectively, we give below the explicit ex- 
pressions for Mi"(p, ki,k 2 )- 

The amplitude corresponding to the non-resonant part a) 
is written, following Novikov and Shifman in JD], as 



Mi" 




(4) 



/3(AM) 2 (l + ^)+ 



|/3((AM) 2 -mL)(l-g)(c OS ^-i)], 

Here AM = M Y „ - M r , = 130 MeV is the pion de- 
cay constant, = \j {k[ + k 2 ) 2 is the invariant mass of 
the two outgoing pions, and 9 is the angle between the vr~ 
and Y b in the dipion rest frame. Eq. © is a guess to model 



the vrvr continuum, inspired by the decay characteristics 
of the dipionic transitions involving Quarkonia states |3|], 
such as T(4S) — > T(lS')vr+vr _ , in which the dipion mass 
spectra do not show any resonant contributions. However, 
as we show here, the dynamical quantities F (a form fac- 
tor) and /3 (a measure of D-wave contribution) required to 
fit the data from the decays Y b — > T(1S, 2S')vr+vr~ are 
very different in magnitude from those required in the de- 
cay T(45) -> T(15)vr+vr- H. 

The amplitude -M^ coming from the diagram b) is 
the resonant part involving the 0++ states / (600) and 
./o (980), and the subsequent decays / (600), / (980) -4 



(5) 



where fo(i) are the two 0++ resonances and the various 
dynamical factors are defined below in terms of the relevant 
vertices and the propagator: 



5/0(0^1-^2 , 




fc2 " m ft(i) + im f0 (*) r fo(i)(' m T"') ' 

(6) 

and /„(*) = / (600) or / (980). The couplings 
5/o(600) = -c/ and g fo{98 o) = V^Cj are taken from |H, 
where c f = 0.02 ± 0.002 MeV' 1 and cj = -0.0025 ± 
0.0012 MeV~ . We use the central values for the cou- 
plings. The propagator of / (600) should not be taken 
in the minimal width approximation, since the total de- 
cay width and the mass are of the same order lllll, [l5fl. 
Following [16], the width is multiplied by a momentum- 
dependent factor: 

«Vo(600) P* 



r(m„ 



/o(600)- 



Po 



(7) 



where = p*{m fo{eoo) ) andp* = p*(m„) are the de- 
cay momenta in the resonance rest frame. The other scalar 
(/ (980)), having r /o(980) /m/ o(980) <C 1, is taken in the 
minimal width approximation, i.e. r(m 7r7r ) = ry (g 80 ). 
The amplitude Ai^ coming from diagram c) is 



M!? = g*'A Ml270) (Tn mr )=g 



a , V8vr(2J+l) y2 



«/ 2 (1270) V m /2 (1270) 



m 



/a (1270) ~~ m ™ ~ ^ m /2(1270)r /2 (i270) 



(8) 



For / 2 (1270), J = 2, and we have kept only the helicity- 
2 component of the D wave with Y 2 2 the corresponding 

spherical harmonics, |Y" 2 2 | = J sin 2 9. In principle, 

there is also a helicity-0 component of the D wave Y° 
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present in the amplitude, but following the high statis- 
tics experimental measurement of the process 77 — > 
/ 2 (1270) -> 7T+7T - by Belle U3], this contribution is 
small, characterized by the value of r 02 , the helicity 0- 
to helicity 2 ratio in / 2 (1270) —> tttt, r 02 = (3.7 ± 
0.3t 2 5 9 9 )%. This can be included as more precise mea- 
surements become available. 

The described diagrams yield a coherent amplitude, and 
the various contributions interfere with each other having 
non-trivial strong (interaction) phases, which are a pri- 
ori unknown. We treat them as free parameters to be 
determined by the fits to the Belle data. Combining all 
three amplitudes, the complete decay amplitudes for Y b — > 
T(15) 7T+7T-, T{2S) vr+vr- are: 



M 




Y T 
£ .£ 



L ![™L-/3(AM) 2 (l + ^) + 



/3((AM) 2 -mL)(l 



2 

+E 



Ami 



mt 



)(cos 2 



o A (o^ fcW K-K)/2 



mi 



m /o(i)+» m /»(i) r /o(i)K) 

+a h (i 27 o) e^^ <i270) A 



(9) 



where <Z/ (j) = gf Q a\Ff Q u)F v , The sum over i runs over 
all ++ resonances contributing in the given energy range. 

The differential decay width (averaged over the polar- 
izations of the initial Yi,-hadron and summed over polar- 
izations of the final T-meson) is given by 

1 1 



dT 



(2vr) 3 32M3 b 



\M\ 2 dm 2 Tn dml„, (10) 



where m^, 



{p+ki) 2 (the amplitude is symmetric under 
the interchange of the two pions). The cos 9 dependence is 
given implicitly by m T7r . By integrating over the phase 
space, we derive the two distributions in m„ and cos 9. 

We have undertaken fits of the Belle data [1] with 
our model (|9), normalizing the distributions for the 
T(1S)tt + tt^ and T(2S)tt + it~ channels to yield the mea- 
sured partial decay widths r T (i S ) +27r = 0.59 db 0.04 ± 
0.09 MeV and r T(2S)+27r = 0.85 ± 0.07 ± 0.16 MeV. 
The input parameters given in Table J] are taken from the 
PDG 111 ill , except for the f (600), for which we have taken 
the values from E791 Ml . 



TABLE I. Input masses and decay widths (in GeV) of the reso- 
nances /o(600), / (980) and / 2 (1270). 



M Yb 


10.890 


m /o(600) 


0.478 


r /o(600) 


0.324 




9.460 


m /o(980) 


0.980 


r /o(980) 


0.07 


M T (2S) 


10.023 


m / 2 (1270) 


1.270 


r / 2 (1270) 


0.185 



The dipion invariant mass distribution dT/dm„ n and 
the angular distribution dT /d cos 9 [GeV] measured by 
Belle [1] for the final state Y(2S)tt + it~ are shown in Fig. 
Q] The shaded histograms are the corresponding theoretical 
distributions from our model having a x 2 /d.o.f. ~ 9/8 
(obtained for the m ww spectrum), with the fit parameters 
given in Table HH yielding an integrated decay width of 
T(Y b -> T^SV+tt-) = 0.85 MeV. The solid curves are 
the distributions for (3 = from the non-resonant part (|4]l 
alone, which are the anticipated distributions from the de- 
cays T(55) -> T(25)vr+7r- JKH. The dashed curves 
correspond to the best- fit solution without the f (600) con- 
tribution, yielding f3 as 0.4 with x 2 /d.o.f w 23/10 (ob- 
tained for the m^Tt spectrum). The difference between the 
histograms (our fits) and the curves is that the latter do not 
have the / (600) contribution. Both the solid and dashed 
curves fail to describe the Belle data. 



0.004 




dr/dcos6 [GeV] 

0.0008 [ 



1 



/ 

k. 



cos 6 

FIG. 1. Dipion invariant mass (m^Tr) distribution (left frame) and 
the cos 9 distribution (right frame) measured by Belle 1 1] for the 
final state r(2S> + 7r~ (crosses), and the theoretical distributions 
based on this work (histograms). The solid and dashed lines show 
purely continuum contributions for different /?. 

TABLE H. Fit values, yielding F = 0.86 ± 0.34, /3 = 0.7 ± 0.3 
for the non-resonant contribution, and for the parameters en- 
tering in the resonant amplitude from /o(600) for the decay 

Y b ->■ T(2S>+7r-. 





a Mi) 


F h(i) 


F Mi)/ F 


tPMi) (rad.) 


/o(600) 


10.89 ± 2.4 


4.19 ±0.92 


4.86 ±2.18 


2.76 ±0.22 



The measured spectra (in m m and cos 9) for the final 
state T(15)7r + 7r _ from Belle [1] are shown in Fig. |2]to- 
gether with our theoretical distributions (histograms) ob- 
tained for the model in (O having a x 2 /d.o.f. ~ 5/5 
(obtained for the ra„ spectrum in the upper left frame), 
with the fit parameters given in Table Hill yielding an in- 
tegrated decay width of T(Y b -)■ T(2S)tt + tt-) = 0.66 
MeV. The two curves in the upper frames show the shape 
of the continuum contribution based on @, with the solid 
curves obtained for f3 = (as would be expected for the 
transition T(5S) — > Y(l,S)7r + 7r~), and the dashed curves 
corresponding to the best-fit solution without the resonant 
contributions yielding (3 ~ 0.3 with x 2 /d.o.f ~ 65/11 
(obtained for the m„ n spectrum). Both of them fail to de- 
scribe the Belle data. In addition we show the contributions 
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from the continuum plus a single resonance in the lower 
frames (solid curves: / (600) with x 2 /d.o.f ~ 16/9; 
dashed curves: / (980) with x 2 /d.o.f « 30/9; dotted 
curves: / 2 (1270) with X /d.o.f « 33/9). They also fail 
to describe the Belle data. 



dr j dm 



dr/ d cos [GeV] 




1.0 

'. GeV ] cos 

FIG. 2. Upper Frames: The distributions measured by Belle (]]] 
for the final state T(lS)7r + 7r _ (crosses), and the theoretical dis- 
tributions based on this work (histograms). The solid and dashed 
lines show purely continuum contributions for different /?. Lower 
Frames: Contributions with continuum plus a single resonance 
(solid curves: /q(600); dashed curves: /o(980); dotted curves: 
/a(1270)). 

TABLE HI. Fit values, yielding F = 0.19 ± 0.03, /3 = 0.54 ± 



0.12 for the non-resonant contribution, a 



/ 2 (1270) 



= 0.5 ± 0.16, 



V/ 2 (i270) = 3-33 ± 0.06 for /2(1270), and for the parameters 
entering in the resonant amplitude from /o (600) and /o (980) for 
the decay Y b -> T(lS> + 7r~. 





a Mi) 


F fod) 




<Pf (i) (rad.) 


/o(600) 


3.6 ±0.7 


1.38 ±0.27 


7.34 ± 1.94 


1.14 ±0.14 


/o(980) 


0.47 ±0.02 


1.02 ±0.04 


5.42 ± 1.0 


4.12 ±0.3 



We also remark that using the fits of the data for the de- 
cay Y b — > Y(1S)tt + it~ presented here, we are able to ex- 
plain the decay width for the decay Y b —> T(1S)K + K~ , 
measured by Belle yj]. The decay is anticipated to be 
strongly dominated by the ++ tetraquark state / (980). 
Details will be published elsewhere. 

Summarizing, we have argued here that the decays Y b — > 
T(1S, 2S)ir + ir~ are radically different than the similar di- 
pion transitions measured in the T(4S) and lower mass 
Quarkonia. The dynamical model presented by us will be 
tested in great detail with improved data, which we expect 
in the near future from Belle. 
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